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Abstract 

The velocity field and the associated tangential stress corresponding to flow of a general- 
ized second grade fluid between two infinite coaxial circular cylinders, are determined by means 
of the Laplace and Hankel transforms. At time t = the fluid is at rest and at t — + cylinders 
suddenly begin to rotate about their common axis with a constant angular acceleration. The so- 
lutions that have been obtained satisfy the governing differential equations and all imposed initial 
and boundary conditions. The similar solutions for a second grade fluid and Newtonian fluid are 
recovered from our general solutions. The influence of the fractional coefficient on the velocity of 
the fluid is also analyzed by graphical illustrations. 

1 Introduction 

A large class of real fluids does not exhibit the linear relationship between stress and the rate 
of strain that is now in great interest of scientists and engineers. Generally, rheological properties 
of a material are specified by their so called constitutive equations. Among the many constitutive 
assumptions that have been employed to study non-Newtonian fluid behavior, one class that has 
gained support from both the experimentalists and the theoreticians is that of Rivlin-Ericksen 
fluids of second grade. The Cauchy stress tensor T for such fluids is given by [1, 2] 

T = -pl + fiAi +aiA 2 + a 2 A?, (1) 

where p is the pressure, I is the unit tensor, fj, is the dynamic viscosity, a\ and ui are the normal 
stress moduli and Ai and A 2 are the kinematic tensors. In the last years, many authors have made 
use of rheological equations with fractional derivatives to describe the properties of fluids. The 

* Corresponding author: Amir Mahmood 
E-mail address: amir4smsgc@gmail.com 

Permanent address : Abdus Salam School of Mathematical Sciences, GC University 
Lahore, PAKISTAN. 



1 



constitutive equations with fractional derivatives have been proved to be a valuable tool to handle 
viscoelastic properties. In general, these equations are derived from known models by substituting 
the time ordinary derivatives of stress and strain by derivatives of fractional order. 

The constitutive equation of the generalized second grade fluids has the same form as (1), but 
the kinematic tensor A 2 is defined by [3-5] 

A 2 = Df Ai + Ax (gradv) + (gradv) T Ai , (2) 

where v is the velocity field, Ai = gradv + (gradv) T , the superscript T denotes the transpose 
operator, and Z?f is the Riemann-Liouville fractional derivative operator defined by [4] 

In the above relation T(-) is the Gamma function. This model reduces to the ordinary second 
grade fluid when = 1, because D\j = df/dt. 

In this paper, we study the motion of a generalized second grade fluid between two infinite 
concentric circular cylinders, both cylinders are rotating around their common axis (r = 0), with 
constant angular accelerations. By means of the Laplace and Hankel transforms we obtain the 
velocity field and the adequate shear stress. 

2 Rotational flow between concentric cylinders 

Let us consider an incompressible second grade fluid at rest in an annular region between two 
straight circular cylinders of radii i?i and R 2 (> Ri)- At time t = + , both cylinders suddenly 
begin to rotate about their common axis with constant angular accelerations. Owing to the shear, 
the fluid is gradually moved and its velocity in cylindrical coordinated (r, 6, z) is given by [2,6] 

v = v(r,t) = w(r,t)ee , (4) 
where eg is the transverse unit vector. The basic equations corresponding to this motion are [6,7] 

r(r,t) = ( M + aiDf)(|:-jMr,t) ) (5) 



^M = (v + aDf)(^ + i|--^)a,(r,t), r € (Ru Jfc), t > 0, (6) 

where r{r,t) = S r $(r,t) is the shear stress which is different of zero, v = p/p is the kinematic 
viscosity, p is the constant density of the fluid and a = ot\j p. The appropriate initial and boundary 
conditions are 

w(r,0)=0, (7) 

u(Ri,t) = JJifiii, uj(R 2 ,t) = R 2 n 2 t fori>0. (8) 
To solve this problem, we shall use as in [7,8] the Laplace and Hankel transforms. 
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2.1 Calculation of the velocity field 

Applying the Laplace transform to Eqs. (6)- (8) and using the Laplace transform formula for 
sequential fractional derivatives [4] , we obtain the following ordinary differential equation 



(u + aq 13 ) 



d 2 u)(r, q) 1 duj(r, q) u(r, q) 



- qU(r, q)=Q, (9) 



dr 2 r dr r 2 
where the image function Z0(r, q) — J °° w(r, t)e~ qt dt of u)(r, t) has to satisfy the conditions 

uJ(R 1 ,q) = ^, uJ(R 2 ,q) = ^, (10) 

q being the transform parameter. 

We denote by ZJ H (r n ,q) — ru(r, q)Bi(rr n )dr, the Hankel transform of the function U(r, q), 
where 

B\{rr n ) = Ji(rr n )Y 1 (R 2 r n ) - Ji(R2r n )Y 1 (rr n ), 

and r n are the positive roots of the transcendental equation Bi(R\r) = and Ji(-) and Yi(-) are 
Bessel functions of order one of the first and second kind. Applying the Hankel transform to Eq. 
(9), taking into account the conditions (10) and using the following relations 

^-[Si(rr n )] = r n [Jo(rr n )F!(ii 2 r n ) - Ji(R 2 r n )Y (rr n )] - -Bi(rr„) (11) 
dr r 

and 

J (z)Y 1 (z) - Ji(z)Y (z) = -—, (12) 

we find that 

/ , 0, f 2 [^2 Jl (glT„) ~ gigi Jl{R 2 r n )] 2 _ \ _ 

L Ttq J\ \R\v n ) J 

or equivalently 

_ . > _ 2[i? 2 n 2 Ji(i?ir n ) - gi^Xggrn)] i/ + agg , . 

^ r "' 9j - ^Ji(i?ir„) q>[q + arlqP + vrlY 

Eq. (13) can be written in the following equivalent form 

(Trio q), (14) 



where 



and 



_ 2 R 2 fl2Ji(Rir n ) - RiCliJi(R 2 r n ) 1 

^ (r "' 9) = ^ Ju^o ? ' (15) 



_ f 2 R 2 fl 2 Ji(Rir n ) - RiVLiJi{R 2 r n ) 1 
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The inverse Hankel transforms of the functions u 1H and ui 2H are 



Wl(r ' 9) = {Ri~Rl)r ' 



respectively, 



-, \ 7T \ - rlJ^R^B^rrn) _ 



Now, we find that the function w(r, q) has the form 

fiii?? - r 2 ) + f! 2 i? 2 (r 2 - i?f) 1 



u(r, q) = 



n=l 

We introduce the notation 



{Rl - R\)r 

^niJi(fl 2 rn)]Bi(rr n ) 
J 2 (i?ir n ) - Jl{R 2 r n ) q[q + ar^ + vr*] 



^ Ji(flir n )[i? 2 n 2 Ji(flir n ) - fliniJiCR 2 rn)]fli(7Tn) 1 , . 

l?/R, r „U./?(J?„rJ «r«-l-™-2„/8 j. ,,,.21 • ^ y ^ 



= F ~ 2V7 21 , (20) 
q[q + arlqP + z^r 2 ] 



and rewrite Eq. (20) in the equivalent form 

,-1-/3 ~ „-l-P-k0 



(q 1 -? +arl) + vrlq-P £ ( ^ {q^ + ar* Y+ 1 ' (21) 



In order to determine the inverse Laplace transform of the function u>(r, t) we will use the following 
formulae [9] 

i ~ 1 (-r} = S r T' a>0 ' 

U a J r(o) 

X 1 (F^ / = Ga - bM t] = § r (c) r ( , + 1) r [(c + ,> - b] ' & )>°- 

So we find that the velocity field u(r, i) has the following form 

_ ^(fl 2 - r 2 ) + r> 2 fl 2 2 (r 2 - gf) 
U ^ t >~ (Rl-R\)r 

^ Ji(.Rir n )[JZ 2 n 2 Ji(fiiT- n ) - fliniJi(fl 2 r n )]Bi(rr n ) 
^ J?(Rir n )-J?(R*r n ) 



x ^(-^ 2 ) fc G 1 _ /3 ,_ 1 _ /3 _^, fe+1 (-ar 2 , t) , (22) 

fe=0 
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or, the equivalently 



^_ ^(fl|-r 2 ) + » 2 ^(r 2 -i^ 
UJ ^ t) ~ (Rl-R\)r 



-E 



n=l 



Ji(flir n )[fl 2 n 2 Ji(flir n ) - RMi{.R 2 r n )\B^rr n ) 
J?(Rir n ) - Jl{R 2 r n ) 



xE 

j,k=0 



{-vr 2 n ) k {-arlyT{k + j + 1) f (i-W+fc+i 

r(fc + i)r(j + 1) r[(i - p)j + k + 2] ' 



(23) 



2.2 Calculation of the shear stress 

The shear stress r(r, t) is obtained from Eqs. (5) and (23) by means of the Laplace transform. 
From Eq. (5) we find that 



d 1 

z { r , q) = (fJ> + aiq (3 ){- )ZU(r, q) . 

or r 



(24) 



Now, applying the Laplace transform to Eq. (23), differentiating the result with respect to r and 
using Eq. (11), we obtain 



2RjRj{Sl 2 -n 1 ) t 1 



^) + 7r E 



n=l 



-Bi(rr n ) - r n B(rr n ) 



x Ji(flir n )[fl 2 n2 Ji (gir„) - gigi Ji(fl 2 r n )] (-^ 2 ) fc (-mj )'T(fc + j + 1) ;; 



J?(iiir n ) - J?(R 2 r n 



E 

j,fe=0 



r(fc + i)T(j + 1) 



i 



,(l-/3)j+fc+2 



+ Ql- 



1 



,(l-/3)j+k+ 2 -/3 



where 



(25) 



B(rr n ) = Jo(rr„)yi(i? 2 r„) - Ji(R 2 r n )Y (rr n ) . 
Applying inverse Laplace transform to the image function r(r,q), we find the shear stress 



r(r,t) = 



{Rl - R*y 



a.\t 



1-/3 



r(2 - p) 



n=l 



-Bi(rr„) - r n B(rr n ) 



Ji(Rir n )[R2^2Ji(Rir n ) ~ RiniJi{R 2 r n )] ^ (-^) fc (-a^)JT(fc + j + 1) 

Jl(R irn ) - Jl(R 2 r n ) .j^ r(fc + i)r(. 7 + i) 



t (l-0)j+k+l 



+ U\ 



t (l-0)j+k+l-0 



r[(i - p)j + k + 2] r[(i - /?).? + fc + 2 - /?] 



(26) 
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3 Limiting case ((3 = 1) 

Making = 1 into Eq. (22) we obtain the velocity field 



-E 



ra=l 



Ji(i?ir n )[i? 2 n 2 Ji(i?ir n ) - fliniJi(fl 2 r n )]ffi(flirn) 
^(i2ir„) - J 2 (i? 2 r„) 



(27) 



corresponding to an ordinary second grade fluid, performing the same motion. Similarly, from 
(26), we obtain the associated shear stress 



n=l 



-Bi(rr„) - r n B(rr n ) 



Jl{R irn )-JKR 2 r n ) 2^ r(fc + i)r(. 7 + i) 



t k+l 



+ ai- 



r(k + 2) r(& + i)_ 

The above relations can be simplified if we use the following relations: 



(28) 



£(-^) fe G ,-2- fe , fe+ i(-m-2,i) = J2(-^n) k 



k=0 



k=0 



^ (-ariyT(k + j + l) t k+l 

r(fc + i)r(j + i) r(fc + 2) 



= E 



oo 

= -— V L 

„ r 2 1 1, i 



(-i/r^) fc t fc+1 1_ 

^ o r(fc + 2) (l + ar 2 )^"" vr l^(k + l)\\ 1 + arl 



1 — exp 



vrit 



1 + ar 2 

As a result, we find that, the velocity field has the form 



_ ^ 2 (^-r 2 ) + » 2 ^(r 2 -^) 



7T y, ^(i?!^)^^^^^) - iZifli Ji(R 2 r n )} Bi(rr„) 



n=l 



J 2 (i?ir„) - J 2 (i? 2 r„) 



1 — exp — 



1 + or 2 



(29) 
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and the shear stress has the form 



r(r, t) 



{Rl - Rfy 



{jit + ai) + 7r 



n=l 



-Bi(rr n ) - r n B(rr n ) 
r 



Ji(Rir n )[Rail2Ji(Rir n ) - RiSliJi(R2r n )] f A* 



Ji(R irn ) - J 2 (R 2 r n ) 



vrt, 



1 — exp — 



1 + ari 



+ 



ai 



1 + ad 



exp 



1 + ar^ 



(30) 



Eqs. (29) and (30) are identical with those obtained by Fetecau et al [6 , Eqs. (3.12) and (3.16) 
for A — > 0]. Making a — > into Eqs. (29) and (30), the similar solutions corresponding to the 
Newtonian fluid, performing the same motion, are recovered. Making Oi = and Q 2 — or 
ill = Jl and O2 = into Eqs. (23) and (26), we obtain the velocity field and the adequate shear 
stress corresponding to the flow between two cylinders, one of them being at rest. 



4 Conclusion and numerical results 

In this paper we establish exact solutions for the velocity field and shear stress corresponding to 
the flow of a generalized second grade fluid between two concentric circular cylinders. The motion 
is produced by the two cylinders which at time t = + begin to rotate around their common axis 
with angular velocities flit and Q 2 t. The solutions, obtained by means of Laplace and Hankel 
transforms, are presented under integral and series forms in terms of the generalized G-function, 
and satisfy all imposed initial and boundary conditions. For (3 = 1 or (3 = 1 and a — 0, the similar 
solutions for the ordinary second grade fluids, respectively, Newtonian fluids are recovered. The 
velocity field and adequate shear stress corresponding to the flow between two cylinders, one of 
them being at rest, are obtained as particular cases of our general solutions. Making ill = and 
0,2 = O into Eqs. (23), for instance, we obtain the velocity field 



QR 2 2 (r 2 - Rl) + do ^ Jl(R 1 r n )B 1 (rr n ) 



(Rl-R 2 )r ^ i J 2 (R 1 r n )-J 2 (R 2 r n ) 



x 



00 



, {-vr 2 n ) k {-ar 2 jY{k+ 3 + 1) fd-W+fc+i 

,j^ o r(fc + i)r(j + i) r[(i-/3)j + fc + 2]' 1 ) 

corresponding to the flow between cylinders, the inner cylinder being at rest. 

Finally, the numerical results are given to illustrate the influence of the fractional parameter 
j3 on the velocity u>(r,t). In all figures we considered R\ = 1, R 2 = 4, Sli = 3, O2 = 1-5, p = 1260, 
ai = 11.34 and fi = 1.48. 

In Figs. 1 the profiles of the velocity u>(r, t), corresponding to the motion of Newtonian fluid 
(the curve wN(r)), second grade fluid (the curve wSG(r))) and generalized second grade fluid (the 
curves wl(r), w2(r) and ui3(r)), are plotted, for different values of the fractional coefficient (3 and 
time t. It is clearly seen from these figures that velocity increases when the fractional coefficient 
decreases. Moreover, the influence of (3 is more strong near boundary of the domain and the 
generalized second grade fluid flows faster than the second grade and Newtonian fluids. 

Figs. 2 depict the histories of the velocity field iv(r,t) at the positions r = 1.3, 2.5 and 3.8, 
for t G [0, 10] and different values of (3. One can see that the influence of (3 is more strong near 
boundary of the domain and the velocity increases when (3 decreases. The units of the parameter 
in Figs. 1-2 are from SI units and the roots r n have been approximated by n7r/(i?2 — Ri) [10]. 
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Figure 1: Velocity profile u(r) for different values of the fractional coefficient (5 
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Figure 2: Time variation of the velocity 
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